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Brief Summary 


"Further Improvements in Program to Calculate Electronic 
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I. LIST OF TASKS ACCOMPLISHED 

I. Characterization of Donors and Acceptors 

The study of donors and acceptors in narrow gap 
semiconductors is by no means simple. Misplaced Hg in MCT seem 
to be the cause, of intrinsic doping. Hg vacancies can produce p 
type materials and Hg interstitials n type. Indium doping has 
caused p type MCT to become n type as expected. It is possible to 
have both deep levels and shallow states. The deep levels can be 
caused by interstitial? and if we neglect them we are left with the 
shallow levels that,|at least for' small x, can merge with the 
conduction band. We can understand this merger as related to the 
Mott transition. When the density of defects increases, discrete 
energy levels become broadened and produce energy levels that 
may join with the energy band. In fact it has been predicted that 
Hg 0 8 Cd 0 2 Te with a donor concentration of at least 10 14 cm -3 would 

not have bound states and this seems to be well confirmed 
experimentally. In general there appears to be no carrier 
freezeout for bandgaps in MCT which have x < .2. The following 
references seem to be relevant. 
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2. Electron-Longitudinal Optic Mode Scattering 

We have studied the Random Element Isodisplacement Model 
of two-mode behavior. We have also reviewed models for 
establishing criteria for two mode versus one mode behavior. We 
have studied effective charges and their meaning for scattering 
electrons. We have also reviewed the literature. All of these 
things are discussed in more detail below. 

A. Szigeti, Transverse and Callen Effective Charge 

Ionic polarizability is caused by the displacements of the ions 
as a whole while electronic polarizability is caused by the 
displacement of the electron cloud relative to the nuclei. The 
displacement of the ions changes the interatomic distance which in 
turn affects the distribution of the electron cloud. This effect is a 
short range interaction, since it only depends on the displacement 
of the near neighbors, There is another type of interaction 
between atomic and electronic polarizations through the internal 
field. This effect is a long range effect since the distant dipoles as 
well as the near dipoles contribute. 

Szigeti defined a correction factor s for deviation from 
heteropolar behavior as a result of mutual distortions of 
neighboring ions due to their overlapping. Szigeti charge represents 
the short range interaction of the electronic and atomic 
displacements. One might expect to find s very different from one 
if the ions penetrate considerably into each other or if the bond 
has partial homopolar character and its polarity changes rapidly 
with the distance between the atoms. Szigeti has shown in his 
paper how the factor s can be obtained in terms of values that 
can all be experimentally measured. 

The transverse charge is defined to give the local polarization 
(electronic and atomic) caused by relative displacements. The 
polarization includes that due to ionic motion and that induced in 
the electronic charges in the absence of electric fields. The 
resulting field which scatters electrons is accounted for by 
separately dividing by e 0 (giving rise to the definition of the Callen 

charge). 


SUMMARY 


e s : The Szigeti Charge e s = se 

This is the effective charge used in figuring the force on an 
ion. Its value is affected by distortion due to overlap. It is a 
measure of deviation from ideal heteropolar behavior. 


e T : The Transverse Effective Charge 



This is defined so as to give local polarization induced by 
relative displacement. We are interested in induced polarization in 
the absence of electric fields. 


e c ; The Callen Effective Charge 


c c = e-r/ea, = 


+ 2 


3e, 


Note for rigid ions e w = 1 and e c = e T = e s . 


B. One and Two Mode Behavior 

In binary alloy systems, lattice vibration spectra can show 
two different behaviors. The first in “one mode behavior” in 
which the phonon frequency changes continuously from the value 
that corresponds to one end member to that of the other end 
member with an approximately constant strength. The second is 
“two mode behavior” in which two frequencies occur with a 
strength approximately equal to the fractional composition of each 
component. 

Experimentally, this behavior can be seen through infrared 
and Raman spectroscopy. A Kramers-Kronig dispersion analysis 
can be applied to reflection spectra to get the real part of the 
dielectric constant Re e(oo), the imaginary part of the dielectric 
constant Im e(oo), and the imaginary part of the reciprocal of the 
dielectric constant -Im (l/etuO). From these functions the 
frequencies of the transverse optical phonons (TO) and those of 
longitudinal optical phonons are determined (LO). The frequencies 
of the TO phonons are determined from the positions of the 
maxima of the Im e(oo) curve and the frequencies of the LO 
phonons from the positions of maxima of -Im (l/e(w)) or from the 
points where the Re e(u>) curve crosses the co - axis [8], In the 
case of the two mode behavior, the Kramer-Kronig analysis yields 
for each x - value two frequencies of the transverse optical 
phonons (T0 1 and T0 2 ) and two frequencies of the longitudinal 

optical phonons (LOj and L0 2 ). 


C. Phenomenological Models Used to Describe Frequency 
Variation in Mixed Crystal Systems: 

i. Virtual crystal: used to treat “one mode behavior” [1,2]. In 
this model all masses and spring constants are taken as 
averages weighted by the mixed crystal composition. 

ii. Cluster model: Verleur and Barker [3], have accounted for 
the two mode behavior in some mixed crystals using a model 
based on short range clustering. 

iii. Random Element Isodisplacement (RED model: the REI model 
which was proposed by Chen [4] and developed by Chang and 
Mitra [9] and also by Harada and Narita [10] who assume that 
in the mixed crystal ABj_ x C x a fraction, 1-x of the immediate 

neighbors of the A-ions are B-ions, and a fraction, x, are C- 
ions. The B ions and C ions always have A-ions as nearest 
neighbors. It also assumed that the atoms of the same kind 
are displaced from equilibrium in the same phase and with 
the same amplitude. Many details about the application of 
this model to compound semiconductors to calculate co(x) is 
given in reference [8], It was shown in reference [11] that 
with the experimental data for (w(x)) the two mode behavior 
in Hgi_ x Cd x Te agrees very well with the theoretical 

calculations using the REI model. In [10] the reflectivity 
experimental data were compared to the theoretical one using 
the REI model. Genzel et al. [12], presented a model 
completely defined by the macroscopic parameters of the 
pure end members. They adopted a type of REI model taking 
into consideration only nearest neighbor interaction. 

iv. Linear chain model [17]: this model was suggested to be the 
best one [5] for providing a criterion for determination in 
advance if a mixed crystal will have one or two mode 
behavior [5]. In this model only nearest neighbor force 
constants are considered. Mazur et al. [6], have discussed the 
kind of impurity modes occuring as a result of isotopic 
substitution in either constituent. This substitution is 
characterized in terms of a term e defined as 


If M* < M, then 1 > € > 0 whereas if M' > M, e < 0. The results 
of the impurity mode calculation are given in ref. 6. In 
substituting for the lighter mass m, a lighter impurity (e > 0) we 
get a localized mode which rises out of the top of the optical band 
and for a heavier impurity (c < 0) we get an optical branch gap 
mode that falls out of the bottom of the same optical branch. In 
the case of substituting for the heavier mass, if we substitute with 
a lighter mass U > 0), we get two modes, a local mode rising out of 
the top of the optical branch and a gap mode rising out of the top 
of the acoustical branch. When substituting with a heavier 
impurity (e < 0) for the heavier mass, no new gap or localized 
modes are generated. According to the study done by Lucovsky et 
al. [5], two -mode behavior would be observed in the crystal 
ABi_ x C x (with m c > m b ) if the substitution of C for B in AB 
produces a gap mode and if the substitution of B for C in AC 
produced a localized mode. For a one-dimensional chain, these 
local and gap modes will occur if A is heavier than B, provided 
that the lighter end member compound AB has a gap in its phonon 
spectrum between the optical and acoustic branches. This criteria 
is applied for MCT and MZT. Although experimentally both MCT 
and MZT show two mode behavior only MZT seems to meet simple 
theoretical criteria for such behavior. 

Some experimental data on far-infrared reflection spectra for 
Hgj_ x Zn x Te are given in reference [13]. Reflectivity data on 
Hgi- x Cd x Te are given on reference [14]. In reference [15] calculated 
values of the optical phonon frequencies as function of composition 
for Hg!_ x Zn x Te and Hg!_ x Cd x Te using REI model are given and 
compared with the experimental values. Also the parameters 
needed for fitting the REI model with the experiment for both 
materials are given in this paper. 
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In the following graph we plot mobility versus 
temperature for several different values of transverse charge 
which have appeared in the literature. The main difference 
appears at intermediate temperature where variations of 10^ or so 
can be seen. See also the accompanying table where numbers and 
references to the values of transverse charge are given. 
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1 
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779743.921 

745762.683 

751620 406 

2 

80.000 

474669.724 

493102.113 
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504562.754 
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475571 316 

3 

100.000 

281514.659 

295403.722 

296080.181 

304674.154 

271577 235 

280418.237 

4 

120.000 

180222.338 

190577.704 
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5 
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128833.736 
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13 
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14 
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29636.356 

30549.166 

26942.153 

28058.859 

15 

340.000 

25048.543 

26373.168 
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27225.181 
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25039.744 

16 

360.000 

22488.594 

23628.672 

23673.498 

24365.650 

21631.696 

22480.530 


AUTHOR 

e T /c (HgTe) 

e T /c(CdTe) 

REFERENCE 

Alex 

2.96 

2.35 S.L. Lehoczky, F.R. 
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Reference -♦Lucovsky et al., Phys. 

Rev. B, Vol. 4, Number 4, 1367-1371, 1971. 



3. Review of Alloy Disorder Scattering 

The alloy Hamiltonian can be written, in an approximation, 
as a uniform virtual crystal Hamiltonian with a residual alloy 
disorder potential. The alloy disordering scattering typically gives a 
term proportional to x(l-x). We have not yet studied this in detail, 
but some relevant references are. 
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3. 


4. Other Scattering Mechanisms 

We have noted that scattering by dislocations may be 
important but we have not yet examined this in any detail. 
Zawadski treat's scattering by dislocations in the following paper. 
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5. 


Further Work on MZT 


We should mention that although MZT has several advantages 
as an IR detector (over MCT) primarily because of its increased Hg 
stability it also suffers from a greater difficulty in growing uniform 
crystals. This is because less Zn than Cd is needed for a given E g 

dEg 

and dx for E ’s of interest is greater for MZT thus small errors 

© 

in x are magnified in energy gap fluctuations. 

Some new papers have appeared. We mention 
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Solutions near q s 0,” J. Phys. France 5£L 3223-3232 (1989). 

2. R. Granger and C.M. Pelletier, “Electron Mobility Evaluation in 
Hgi_ x Cd x Te and Hg!_ x Zn x Te with Two Mode Optical Dispersion,” 
Presented at 5th Int. Conf. on II-VI Compounds, Yokagama, 
Japan. 

3. R. Granger, A. Lasbley, S, Rolland, C.M. Pelletier, and R. 
Triboulet, “Carrier Concentration and Transport in Hgj_ x Zn x Te 
for x near 0.15,” J. Crystal Growth 682-688 (1988). 

4 . S. Rolland, A. Lasbley, A. Seyni, R. Granger, and R. Triboulet, 
"Electrical Characterization of as-grown, annealed and 
indium-doped Hgj_ x Zn x Te for x near 0.15,” Revue Phys. Appl. 

24, 795-802 (1989). 

5. E.A. Patten, M.H. Kalisher, G.R. Chapman, J.M. Fulton, C.Y. 
Huang, P.R. Norton, M. Ray and S. Sen, “HgZnTe for very long 
wavelength infrared applications,” J. Vac. Sci. Technol. B9(3), 
1746-1751 (1991). 


6. Experimental Data 

Good comparisons of our calculations with experimental data 
have been done in the Ph.D. thesis of Dr. Wafaa A. Gobba. The 
comparison will be found on the following graph. The abstract of 
the paper will be found in Section III. Our calculation shows better 
agreement with experiment than does the calculation of the 
experimenters. • ■ See the following graph. 


Electron mobility (cm2 /V-s) 


1000000 



Temperature(K) 


Fig. 7. A comparison between theoretical 
and experimental results for MZT electron 
mobility. 



7. Superlattices 


A superlatticc (SL) which has alternating layers of CdTe and 
HgTe has been proposed by Schulman and McGill as a new 1R 
material which can be made by molecular beam epitaxy. The 
first successful molecular beam epitaxy growth of a Hg-Te 
superlatttice occurred in 1982. The band gap depends on the width 
of the quantum wells as well as the barrier width and height. 
Thus the formation of HgTe/CdTe superlattices allow a good control 
of the bandgap which can be a variety from zero to about 1.6 eV. 
We have also done some rough modeling of high mobility devices 
which are based on doping in the barrier to furnish electrons to 
the quantum wells. Some results are given below. Representative 
references are listed. 

There are several comments one could make about the utility 
of superlattices as IR detectors. In some ways superlattices are an 
advantage. For one, the energy gap versus SL well layer thickness 
d w can be better controlled in SL’s than E g versus x in bulk 
crystals - at least for E g ’s of interest. This is particularly true 
because the cut off wavelength for E g ’s of interest is a much more 
rapidly varying function of x for bulk crystals than of the layer 
thickness for superlattices. 

There exists a neat trick in SL’s that can be used to increase 
the number of electrons without reducing the low temperature 
mobility due to ionized impurity scattering. This is discussed later . 
In calculations dealing with electron mobility in superlattices, it 
has been found that the temperature dependence of the mobility is 
greatly affected by interface roughness which leads to big 
complications in making accurate calculations. 

There are other advantages to the HgTe-CdTe superlattices as 
given by Smith, McGill and Schulman. The super lattice tunneling 
currents are reduced compared to bulk and diffusion currents in 
photovoltaic devices made from SL’s are reduced as is band to 
band tunneling currents across the junction. 


HgTe-ZnTe SL’s have also been considered but because of a 6.5$ 
lattice mismatch these become strained layer superlattices. The 
straining can cause additional effects which may be tailored to 
advantage. 
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In the following graphs we have done a calculation which models 
very crudely the way that mobility can be enhanced by using 
layered structures consisting of alternate layers of quantum wells 
and barriers. If one adds donors in the barriers the resulting 
electrons will drift to the quantum well resulting in more electrons 
without more scattering centers, The actual calculation we have 
done keeps the number of electrons constant and reduces the 
number of scattering centers. Thus the low temperature mobility 
(parallel to the layers) would be increased. In the following graphs 
sit = 10 means we have N d electrons at low -temperatures with 
only N d /10 ionized impurities/ volume. Further increase of sit to 
100 produces minimal low temperature gain in the mobility. There 
is no change in the high temperature mobility reflecting the fact 
that ionized impurity scattering is not important in this regime. 
The fact that the effect seems to “saturate” for high electron 
concentration at low temperature would mean that ii impurity 
scattering has become negligible due to decrease in ii scattering 
centers. 
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I hcorellcal Calculations of Electr on Mobility In Narr ow Gap 
Slmii I conductors * 

WAFAA A. GOBBA, Department of Physics and Space Sciences, Florida 
Institute of Technology, Melbourne, FL 32901 

The prototype narrow gap semiconductor we deal with Is Mercury 
Cadmium Tellurlde (MCT) because of Its applicability to making Infrared 
( IR) detectors. Interest has focused on Hg,_ x Cd x Te-because Its band gap Is 
tunable with x, and because It can operate efficiently as an IR detector at 
77K In the atmospheric window of wavelength 8-12 microns. We also do 
calculations on Mercury Zinc Tellurlde (MZT) and find for equivalent energy 
gaps that Its mobility, as a function of temperature, Is similar to ITCT. 
Because of Mercury stability MZT may turn out to be a better IR detector. 
Our calculations are state of the art and Involve the Boltzmann equation, 
Its solution by the variational principle and they Include all relevant 
Interactions The results compare well with experiments. Comments will 
be made about how the author came to specialize In this area and the 
opportunities and challenges that exist for women. 

* Supported by NASA Contract Number NAG-70 1 . 
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Suggested title of session - Abstract Submitted Physics & Astronomy 

In which paper should be for the Indianpolls Meeting or Classification Scheme 

placed: II -VI Semiconductors the American Physical Society Number 75.20 

16-20 March 199 


Electron Mobility In Mercury Zinc Tellurlde Alloys . " WAFA A A. 
GOBBA, J. D. PATTERSON, Dept, of P hysic s and Spa c e Scienc es. F.I.T..and 
51 tEHOCZKY, Marshall Space Flight Cent er - We have calculated the 
mobility of electrons In n-type Mercury Zinc Tellurlde (MZT) and compared 
It to a calculation of the mobility of electrons In n-type Mercury Cadmium 
Tellurlde U1CT ) with nearly the same energy gap and with the same number 
of donors and acceptors. We find for equivalent energy gaps that the 
mobilities In the two compounds (MCT, MZT) are nearly the same. The 
calculations for both MCT and MZT are based on the best set of materials' 
parameters that we can compile from the literature, and the MZT compilation 
Is much more complete than was previously available. Since MZT Is harder 
and structurally more stable with respect to Hg retention than MCT, the 
possibility of equivalent mobility for MCT and MZT Is significant. This 
calculation Is one of the first extensive calculations or the mobility of MZT, 
and we compare It to another by Rolland et al. (Revue Phys. Appl. 24 . 795 - 
802. (1989 )_whlch appears to be less extensive. Our calculation yields very 
good results ns compared to experiment except al very low temperatures 
where Its accuracy Is limited by a lack of knowledge of Ionized Impurity 
concentration. Mobilities In the range of 10^ cm^/vs at 100 K are typical. 
Our calculation Involves scattering of the electrons by longitudinal optic 
phonons, acoustic phonons, Ionized Impurities, holes, and compositional 
disorder. Since not all of these Inter actions can be approximated by elastic 
scattering, the corresponding Boltzmann equation Is solved by a variational 
principle. * Supported by NASA Contract No. HAG-781 
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Enhanc ed Scr eenin g In Dope d Sem ic onductors . 

J.D. PATTERSON and WAFAA A. GOBBA, Dept, of P hysics a nd Sp ace 
Sci en ces. Florida Institute of Technology — Semiconductors with 
shallow donor centers are considered. The donor energy levels are 
assumed to be discrete but screened by electrons from the conduction 
bond. High doping levels are assumed. The number (n)of electrons In 
the conduction band Is Increased by screening which reduces their binding 
energy causing the promotion of more electrons to the conduction band 
and further reducing the binding energy etc. The equations describing this 
effect are analogous to mean field equations. The effect Is temperature 
dependent and near a certain temperature there can be a rapid Increase 
of n with temperature. The abrupt Increase in n is, of course, reminiscent 
of oilier abrupt changes related to screening such as the Moll Transition. 
We predict It should be observable for reasonable values of parameters 

such as In GaAs at T — 20K with 5 x 10 ^ donors/cm^ and for donor 
binding energy of 12 meV. Some other aspects of screening In 
semiconductors will also be discussed. 
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ABSTRACT 


The electron mobility, concentration of electrons, light holes and 
heavy holes, and Fermi energy in n-type Hgi- x Zn x Te (MZT) have been 
calculated and compared with those of Hgi- x Cd x Te (MCT) for the same 
energy gap (but different x), and the same donor and acceptor 
concentration. The results for MZT were found to be very close to those of 
MCT. The calculation relied on solving the Boltzman transport equation 
using variational principles. The processes that were of 
significance included the scattering of electrons by : ionized impurities, 
holes, compositional disorder, acoustical phonons and optical phonons. 

In the process of calculating the results for MZT a table for the latest 
values of all the needed material parameters was obtained by a literature 
search. 

MZT has been considered as a substitute for MCT for possible use in 
infrared detectors due to its relative stability and hardness. Our theoretical 
results for MZT were compared to the experimental results published by 
Rolland et al. and Granger et al. and very good agreement over most of 
the temperature range was obtained. 


The effect of making changes on some of the relevant parameters in 
the calculation was checked. Also, adding neutral impurities to the 
scattering processes and using the most recent published band 
parameters of MCT (rather than the slightly older values) on the electron 
calculation did not produce significant changes. 


The speed of a very involved and long program has been increased 
by changing its code from Basic for an HP 9845 to FORTRAN for a Vax 
11/780 and later a Sun workstation 386i. This program was given to us 
by Dr. S. L. Lehoczky and it was written to calculate the electron mobility, 
the free carrier concentrations, and the Fermi energy for MCT. All of these 
were calculated as a function of temperature, composition, donor and 
acceptor concentration. The program was checked, corrected and fine 
tuned by including the latest values of material parameters for MCT, 
another version was made to run for MZT and neutral defect scattering 
was added to the program. 
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ELEC I NON MONILI IT IN n-lYPE 
lUJi^CdjJe AND l lg 1 . JI Zii )t le ALLOYS* 


J.D. Patterson and Wafa a. A. Qobba 
Department of Physics and Space Sciences 
Flor ida Institute of Teclmolocjy 
Melbourne, FL 32901-0988 U.S.A. 

and 


31. Lehoczl'y 

E375, Space Science Laboratory 
Marshall Space Flight Center, AL 35812 U.S.A. 


ABSTRACT 


V/e have calculated the mobility or elections In n-lype Mercury Cadmium 
Telluride (MOT) and compared it to a calculation of the mobility of 
elections In n-type Mercury Zinc Tellui Ide (MZT) with nearly the same 
energy gap and with the same number of donors and acceptors. We also 
compared the results of the MZT calculation with experiment. We found 
for equivalent enetgy gaps that the mobilities In the two compounds 
(MOT, MZT) were nearly the same. 1 he calculations for botli MCT and 
MZT were based on the best set of material parameters that we could 
compile from the liteiature. Using these parameters, the comparison with 
experiment for MZT yielded very good results. Since MZ I is harder and 
structurally more stable with i espect to I lg t ete title n than MC r, the 
possibility of equivalent mobility for Me r and MZ I is significant. 

This calculation is one of the fir st extensive calculations of the mobility of 
MZT. and we compared it to another which appeal ed to be less 
extensive. Our calculation involves scattering of the elections by 
longitudinal optic phonons, acoustic phonons. Ionized impurities, holes, 
and compositional disorder. Since not all of these Intel actions can be 
approximated by elastic scattering, the con espo tiding Boltzmann 
equation was solved by a var iatlonal principle. We also discuss 


dilections for future woik. 



